Site percolation on planar <£> 3 random graphs 
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In this paper, site percolation on random <3? 3 planar graphs is studied by Monte-Carlo numerical 
techniques. The method consists in randomly removing a fraction q = 1 — p of vertices from 
graphs generated by Monte-Carlo simulations, where p is the occupation probability. The resulting 
graphs are made of clusters of occupied sites. By measuring several properties of their distribution, 
it is shown that percolation occurs for an occupation probability above a percolation threshold 
p c =0. 7360(5). Moreover, critical exponents are compatible with those analytically known for bond 
percolation. 
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I. INTRODUCTION 

Percolation on regular lattices has been extensively 
studied, both analytically and numerically It is now 
firmly established that critical exponents, governing scal- 
ing laws near percolation threshold, depend only on the 
dimension d of the lattice. On the contrary, percolation 
thresholds depend on the precise structure of the lattice. 
The complete knowledge of site and bond percolation 
thresholds for all regular lattices with given dimension 
is a very interesting challenge, from both theoretical and 
experimental points of view. It is also known that there 
exists an upper critical dimension d c — 6. This means 
that for dimension d > d C) all regular lattices belong to 
the same universality class as a family of regular lattices 
containing no loops, called Bethe lattices. This is due 
to the fact that the proportion of closed loops in a large 
regular lattice decreases with the dimension d and even- 
tually becomes negligeable for large d. 

On the other hand, percolation on random graphs is 
still an open subject. One class of such random graphs, 
called complex networks, containing scale free networks 
and Erdos-Renyi networks, is at present attracting a lot 
of interest in physical and mathematical communities as 
they are good models for real networks (world wide web, 
social networks, . . .) [2(], A crucial feature of complex 
networks is that closed loops can be neglected for large 
graphs. Percolation theory have been recently used to 
investigate their intrinsic properties 0|. 

A radically different family of random graphs have 
been extensively studied in the past decades as a non- 
perturbative regularization of quantum gravity (see [3| 
for a review). Contrary to complex networks, they are 
planar, with closed loops that cannot be neglected. More- 
over, distant vertices are strongly correlated. In some 
sense, these graphs are closer to regular lattices than 
complex networks. More precisely, they look locally like a 
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regular lattice but are globally very different. This inter- 
mediate situation makes the problem of percolation on 
these graphs very exciting. Planar $ 3 random graphs, 
or their dual planar dynamical triangulations, belong to 
this family. They are defined in section II of this paper. 
Their properties are rather well established now on the 
ground of a great amount of analytical and numerical re- 
sults In particular, it is known that the Hausdorff 
dimension of these graphs is da =4 0| with a fractal 
structure of so-called baby universes [2J. Percolation on 
planar <I> 3 random graphs is the subject of this paper. 
In fact, bond percolation on planar <& 3 random graphs 
has been exactly solved as the limit q — > 1 of a q-state 
Potts model, using matrix models [§]. In this article, we 
study numerically site percolation on these graphs. The 
main purpose is to measure the value of the percolation 
threshold. Moreover, our work is a test of universality 
between site and bond percolation for this model. 



II. THE MODEL 
A. Random <E> 3 planar graphs 

We consider the set of all planar graphs with N triva- 
lent vertices , i.e graphs without boundaries that can be 
drawn on a sphere and where each vertex is linked to ex- 
actly three neighbors. Morover, two distinct vertices can 
be linked by at most one link and no vertex can be linked 
to itself (see Fig. [Q). 

These graphs are purely topological objects as no 
length scale is given here. Such graphs are characterized 
by their Euler number x = N—Ni+Nf = 2, where N, Ni, 
and Nf are respectively the number of vertices, links and 
faces. Moreover, local properties of these graphs imply 
2N[ = 3N. However, there is no constraint on the size of 
a face, i.e the number of links surrounding a face, except 
that it must be greater than three. Note that, in some 
cases (degenerate graphs), some loops contain almost all 
links of the graph. This set, called $ 3 planar graphs, is 
turned into a statistical model by assigning a Boltzmann 
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Figure 1: A $ 3 graph. 



weight to each graph. In this paper, each graph has the 
same weight. The partition function of this ensemble of 
random graphs is written 



Z N = £ 



1 



Ge* 3 



G(G) 



where the sum is over all <I> 3 graphs with N vertices as 
defined above, and G(G) is a symmetry factor wich avoid 
the overcounting of some symmetric graphs (C (G) is al- 
most always equal to one for large graphs). The size of 
each face can be seen as a random variable, with mean 
value equal to six, whose distribution can be exactly cal- 
culated However, toplogical constraints (planarity) 
imply strong correlations between distant faces. As an 
example, correlation between ajacent faces follow a mod- 
ified Aboav's law [lol |. Another important feature of 
these graphs is their Hausdorff dimension, which has been 
shown to be dn = 4 It can be defined as follows: first 
define a path between two vertices v i and V2 as a succes- 
sion of adjacent links connecting v\ and V2 ■ The length 
of the path is the number of links in the path and the 
geodesic distance between V\ and V2 is the length of the 
shortest path between them. (N r ) a , the mean number of 
vertices whose geodesic distance from an arbitrary vertex 
v a is lower than r, scales (for large r) as 

(N r ) o ~ r d » 

B. Site percolation 

1. Definition 

We now consider the problem of site percolation on 
these graphs. As usual, each site (vertex) of a graph 
G is randomly occupied or empty, independently of the 
rest of the sites. More precisely, each site is occupied 
with probability p or empty with probability q = 1 — p 
( see Fig. [2]). Each distribution of occupied and empty 



Figure 2: Occupied (black bullets) and empty sites on a <3? 3 
graph. Dashed links connect empty sites to occupied or empty 
sites. 



sites on G induces a structure of clusters. A cluster is 
a set of occupied sites connected by links of G. The 
study of average properties of these clusters as p is varied 
is the subject of percolation theory. Consider a graph 
G G $ 3 |jv and a distribution of occupied/empty sites 
on G, denoted V (for a given occupation probability p); 
suppose that a quantity A = A [G, V] depends on G and 
V. The annealed average quantity (^4) (N,p) is obtained 
by first averaging on all distributions V for a given graph 
G and, then, by averaging on all graphs G. It can be 
written: 



(A) (N,p) 



E 

G6* 3 I 



1 



G(G) 



(A) (N 7 p) can be seen as the (average) value of A on a 
graph picked at random from the *1> 3 1 jv ensemble . 



2. Cluster distribution 

For fixed p, the distribution of the sizes of clusters 
is described by the quantity n(s,p), equal to the den- 
sity of clusters made of s connected sites. Alterna- 
tively, one is interested by all the moments of the dis- 
tribution: K n (p) = J2 S n ( s iP) s ™) where the symbol 

means that only finite clusters enter the sums. In 
practice, only the first moments are usually studied. 

More precisely, K (p) = J2 S n ( s >P) ls the total num- 
ber per site of finite clusters; K\{p) = ^ s n(s,p) s is the 
probability that a site belongs to any finite cluster and 

S(p) = ^ s n(s , p) s 2 1 Y^, s n ( s iP) s is ° ne wa Y to define 
the mean cluster size. 

Other interesting quantities are s 1 (p), s 2 (p), . . ., de- 
fined as the sizes of the largest, second largest,. . . cluster. 
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Percolation threshold 



with an exponent r as 



One of the most basic questions is the existence, in the 
thermodynamic limit N — > oo, of a percolation threshold 
p c such that, for p < p c , all clusters have a finite size and, 
for p > p c , there exists at least one spanning cluster. This 
is formally expressed by defining a probability R{p) that 
there exists a spanning cluster, with R(p) = for p < p c 
and R(p) = 1 for p > p c . The definition of a spanning 
cluster is not unique but, practically, it can be seen as a 
cluster connecting two opposite boundaries on a lattice. 
Unfortunately, such a definition does not easily apply for 
the family of graphs considered in this paper, as there 
is no natural way to define opposite boundaries in these 
graphs. However, this can be circumvented by defining 
p c as the value of p for which the sizes of the second, 
third, ... largest clusters reach a maximum [111]. This 
property can be understood with the following heuristic 
argument: for p < p c , all graphs contain several finite 
clusters with comparable size, as soon asp is not too close 
to p c ■ As p is growing, the sizes of the largest, second 
largest, third largest . . .clusters are growing until p c is 
reached. At this point, most of the clusters merge into 
a giant component so that, for p > p c , the sizes of the 
second, third,. . . largest clusters dramatically fall down. 



4- Critical exponents 

There is a profound analogy between percolation and 
critical phenomena, where p c plays the role of the critical 
temperature. The analog of the order parameter is the 
density of the spanning cluster, denoted P(p) . It repre- 
sents the probability for a site to belong to the spanning 

cluster. It can be defined as V(p) = p~Yls n ( s >P) s smce 
p is the probability for a site to belong to a cluster (finite 

or not) and J2 S n ( s >P) s ls the probability for a site to 
belong to a finite cluster. For p < p c , all clusters are fi- 
nite and V(p) — 0. For p > p c ,V(p) increases with p and 
< V(p) < 1. However, for planar <& 3 random graphs, 
the notion of spanning cluster is not appropriate, and 
V{p) is defined as the density of the largest cluster. At 
least for p ~ p c and large graphs, both definitions are 
expected to coincide. 

As in critical phenomena, it is possible to define a cor- 
relation length £ (p) related to the mean radius of clusters 
and diverging at p = p c as ~ (p — p c ) . Moreover, 
most observables are expected to follow scaling laws in 
the region p ~ p c . We define standard exponents a, (3 
and 7 associated with K (p), V(p), S(p) in the critical 
region by 



K [p) ~ (p- 
V(p) ~ (p- 
S(p) ~ (p- 



Pc 
Pc 
Pc 



n(s) ~ s T 

These exponents are expected to follow the scaling re- 
lations 



2 - a = ~f + 2(3=dv 
T = 2 + J~p 



(1) 



Here, d is a parameter characteristic of the model. For 
regular Z)-dimensional lattices with D < 6, d = D is 
the standard dimension of the lattice. For D > 6, i.e. 
above the upper critical dimension, d = 6. For planar $ 3 
random graphs, d = d H = 4, the Hausdorff dimension 
defined above. 

The sizes of the larg est, second largest, . . . clusters 
follow a scaling law [llj 



si 

S2 



(P - Pc) 
- Pc) 



P-vd 



P-vd 



5. Theoretical values of critical exponents for bond 
percolation 

Percolation could be defined by an occupation proba- 
bility on the links of the planar $ 3 random graphs. This 
would give rise to bond percolation. As in critical phe- 
nomena, critical exponents are expected to be universal 
whereas p c should depend on the details of the model. 
So, critical exponents should be the same in bond and 
site percolation. Bond percolation on planar $ 3 random 
graphs has been exactly solved using a random matrix 
model formulation [s| , so that exponents v, a, (3, 7 and r 
are exactly known in this case 

1 15 

v = 1 , a = -2 , /3 = - , 7 = 3, r = — 



expo: 

m 



and, at p = p c , n(s) = n(s,p c ) is expected to scale 



If universality holds, critical exponents for site percola- 
tion should also be given by these values. 



III. THE NUMERICAL EXPERIMENT 



The method 



Nodes removal 



The first scoincidetage is to generate graphs in <E> 3 1 jv 
. We start from a tetrahedron. Then, one face (trian- 
gle) is randomly chosen, a vertex is added inside this 
face and linked to the three vertices of the triangle. This 
procedure is repeated until a polyhedron with N faces 
(triangles) is obtained. This polyhedron is transformed 
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into a $ 3 graph by duality, i.e. each face (triangle) of 
the polyhedron is replaced by a vertex linked to three 
vertices associated to the three adjacent faces of the tri- 
angle. The graph G^°' thus obtained contains N trivalent 
vertices and has the topology of a sphere. Then, starting 
from G^°\ all graphs in <& 3 \n can be generated by using 
standard flips of links called T 2 moves 

Once a graph G G $ 3 |jv has been obtained, one distri- 
bution V(G,p) of occupied and empty sites is randomly 
generated. This is achieved by randomly removing qN 
vertices from G, with q = 1 — p. The pN remaining 
vertices are defined as occupied vertices on G. Then, in- 
teresting quantities are measured. Several distributions 
are generated this way to obtain average properties on 
G for fixed p. Then, another graph is generated by a 
serie of T 2 moves and the removal procedure is applied 
to this new graph, etc. Annealed averages of observables 
are thus computed. 



2. Cluster construction and measured quantities 

For each distribution V(G 7 p) on a given graph G, all 
clusters are constructed using a breadth-first search al- 
gorithm analog to Wolff algorithm . At step n — 1 
of the algorithm, suppose n — 1 clusters ci, C2, . . . , c n _i 
have already been constructed. All corresponding occu- 
pied sites are labeled "visited". Step n first consists in 
choosing an occupied site not yet visited. This site v 
is the root of the cluster c n and is now labeled "visited". 
It is put in a (empty) list Q. The following procedure 
is now applied to Q : for each site v in Q, all occupied 
and not yet visited neighbors of v on G are added to c„, 
labeled "visited" and put in Q whereas v is removed from 
Q. This procedure is repeated until Q becomes empty. 
The cluster c n is thus completely constructed. The algo- 
rithm stops when all occupied vertices have been visited. 

The size of each cluster - i.e. the number of sites in 
the cluster - is registred and a histogram of the sizes is 
built. This allows to measure the following quantities: 

• K (p) , the total number of clusters (except the 
largest one); 

• n(s,p) , the density of clusters of size s; 

■ V(p) , the size of the largest cluster divided by N; 
for large N, this is expected to represent the order 
parameter defined above; 

• S2{p) , the size of the second largest cluster ; 

• K2(p) , the second moment of n(s,p); in the crit- 
ical region, K 2 {p) is expected to scale like S(p) = 
K 2 (p)/K 1 (p), the mean (finite) cluster size . 



3. Finite size scaling 

An important tool when using numerical simulations 
is the finite size scaling analysis. It is based on the hy- 
pothesis that, for finite systems, scaling laws are cor- 
rected by scaling functions depending on the ratio be- 
tween the linear size L of the system and the corre- 
lation length £ . More precisely, for L ^> £, the sys- 
tem should not feel finite size effects, but when the cor- 
relation length becomes of order L, finite size should 
modify scaling laws. For a quantity O, a scaling law 
O ~ (p — Pc) z (for infinite size of the system) can be 
rewritten O ~ S, z ^ as the correlation length itself scales 
as £ ~ (p-pc)~ v - Then, when f ~ L, O ~ L z ' v . This 
can be summarized by O ~ L z ' v F ((p — p c ) L 1 ^) where 
F(x) is a (scaling) fonction of the dimensionless ratio 

x = (p-pc)L x l v ~ such that F(x) -» 1 for 

x ~ 1 and F(x) — ► x~ z for x 3> 1. 

Finite size scaling laws enable one to extract the values 
of critical exponents by studying the behavior of quanti- 
ties as the size of the system is varied. This also defines a 
finite size percolation threshold: suppose that the scaling 
function for O reaches a maximum for x — x Q . Then, for 
fixed L, the value of p for which O reaches a maximum 
is given by (p — p c ) L x l v — x a . This defines an effec- 
tive finite size percolation threshold p c (L) approaching 
p c when L — > 00 as p c (L) — p c ~ L" 1 /" [lj. 

For percolation on graphs in $ 3 |at , there is no explicit 
linear size as graphs are purely topological. However, 
the quantity N 1 ^" , where dn is the Hausdorff dimen- 
sion, defines an effective linear size, so that finite size 
scaling can be written O ~ L z l v F ((p - p c ) N 1 /" 4 ") . 
The effective finite size percolation threshold p c (N) then 
approaches p c as p c (A r ) — p c ~ N~ x l vdH . 



B. Simulations 



We simulated graphs of sizes ranging from N = 400 
to N — 25600 or N = 51200 (according to the measured 
quantity) vertices. When measuring V(p) and S2(p), we 
considered various values of p for each size. For each 
graph G and given p , we generated nx> = 128 to 1024 
occupation distributions by the nodes removal method. 
Tip was chosen big enough to minimize its influence on 
annealed averages. Each simulation consisted in gener- 
ating no graphs in $ 3 |jv, with no ranging from 1024 to 
4096 for each value of p. Each graph was obtained from 
the previous one by performing 2000 Ni flips (T2 moves) . 
We estimated error bars by standard jackknife method. 
Error bars are always plotted on the figures below but, 
most of the time, they are hidden by the symbols. 
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Figure 3: Order parameter V vsp. 



Figure 4: Size of the second largest cluster 82 vs p. 



C. The results 

1. Order parameter 

The percolation order parameter V plotted as a func- 
tion of p is shown in Fig. [3j for various values of N. We 
see that, for low values of occupation probability p, V 
is close to zero. As expected, V = 1 for p = 1 since all 
graphs are connected when no vertex is removed. Ac- 
cording to percolation theory, V is not zero only in the 
percolating phase. As can be seen on the figure, in our 
model, a percolation transition is expected to take place 
for p ~ 0.7. However, it is not possible to get a precise 
estimate of the threshold value p c with these data. 



2. Percolation threshold p c 

The percolation threshold p c is determined using the 
behavior of S2, the size of the second largest cluster. Fig. 
H] clearly shows a peak of S2 for a value depending on 
N, denoted p c {N) . This effective percolation thresh- 
old p c {N) was extracted by fitting the peaks of S2 with 
quadratic functions. The result is plotted in Fig. [H As 
explained above, by using finite size scaling hypothesis, 
p c (N) is expected to approach p c for large N according 
to 



Pc {N) = Pc + Cl N- 



-l/vd f 



(2) 



where c\ is a constant. So, we fitted p c {N) with the 
scaling law ^ and we obtained p c — 0.7360(5). 



3. Critical exponents 
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Figure 5: Finite size percolation threshold vs N. 



The fit of p c {N) with the scaling law ^ allowed us to 
extract also the value of 1/ vdu ■ We obtained 1 / vdn — 
0.489(9). As can be seen in Fig. [6l p c {N) is clearly a 
straight line when plotted as a function of TV -0 - 489 . 

• Exponent j3. 

We also used S2 to extract (i/vdn- Fig. [4] shows that 
the peaks of S2 are growing with the size N. By fitting 
the peaks of S2 with quadratic functions, we obtained the 
value of S2 at the maximum, denoted at s 2 raaa; (iV) . Fig. 
[7] shows the results. 

By finite size scaling arguments, S2 is expected to be- 
have as 
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Figure 6: Finite size percolation threshold plotted as a func- Figure 8: Number of clusters of size s (not normalized) for 
tion of JV -0 - 489 . p = Pc an d N = 51200. The straight line is the best fit. 
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Figure 7: Maximum of S2 us N. The straight line is the best 
fit. 

so that s 2 nax {N) scales as 

s m«i ^ jyl—0/vdH ^ 

By fitting s™ aa: (./V) with the scaling form ([3]), we ob- 
tained 1 - (3/vdn = 0.880(1). The best fit is plotted in 

Fig. m 

■ exponent r 

We measured the number of clusters of size s at p — p c — 
0.7360. In order to extract the value of t, we did not use 
finite size scaling analysis. Instead, we used the scaling 
law (for s>l) 

n(s) ~ s- T (4) 



for one value of N large enough to minimize finite size 
effects. We chose N = 51200 and discarded data with 
s > 100 and s < 7000 (because the scaling form is valid 
for s 3> 1 and there was not enough statistic beyond 
s = 7000). Then, we fitted the data with the scaling 
form Q. We obtained r = 2.0(1). The result is shown 
in FigH 

• Exponent a 

We extracted the value of a by using finite size scaling 
analysis applied to the density of total number of clusters 
K Q (p) measured at p = p c (Fig. [9j. The expected scaling 
law is 

K (N, Pc ) ~ N( a - 2 V udB (5) 

We fitted data with ([5]) and obtained (a — 2) jvdn = 
-0.7(2). The result is plotted in Fig. [9] 

• Exponent 7 

We measured K2 (p) , the second moment of the distribu- 
tion n(s,p) for p = p c . Fig. [10] shows the result. As 
mentioned above, K2(p) is expected to scale as S(p), the 
mean size of finite clusters, at least in the critical region. 
Finite size scaling law for Ki{p) is then 

K 2 (N,p c ) ~ Wl vd » (6) 

We extracted 7 by fitting data with J6]) and we ob- 
tained j/i/dji — 0.67(3) 

4- Remark on the results 

The determination of the percolation threshold using 
the maximum of the second largest cluster gives rather 



7 



0.0580 



0.0578 



o 



0.0576 



0.0574 



1 



10000 20000 30000 40000 50000 60000 

N 



Figure 9: Density of total number of cluster measured at p ■ 
p c - The dashed curve is the the best fit. 
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Figure 10: Second moment of n(s,p) for p = p c . The dashed 
line is the best fit. 



accurate value for p c . This is partly due to the fact that 
the parameter p c in the scaling law (J2j) is independent of 
the precise approach of p c (N) to its infinite size limit. 

The situation is rather different for the determination 
of critical exponents. In Table ITTl measured and theoret- 
ical values of critical exponents considered in this paper 
are summarized. In addition, combinations entering scal- 
ing relations (JXJ are given in Table Hill As can be seen, 
the value of l/vdn obtained by simulations is not com- 
patible with the theoretical value 1/4. In fact, it is well 
known that the Hausdorff dimension of planar <£> 3 random 
graphs is very sensitive to finite size effects and can be 
extracted only for large lattices using sophisticated scal- 
ing variables in finite size scaling analysis i |. However, 
this does not mean that simulations are unable to take 
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(a - 2) /vd H 
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Simulation 


0.489(9) 


-0.7(2) 


0.120(1) 


0.67(3) 


2.0(1) 
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Table II: Theoretical and measured values of critical expo- 
nents. 
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Simulation 
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2.136(1) 
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Table III: Test of scaling relations between critical exponents 



account of the fractal structure of these graphs. The re- 
maining critical exponents are compatible {f3/vdn,T) or 
marginally compatible ( (a — 2) /vdu , l/vdn) with the 
theoretical values. However, error bars take into account 
neither the error on 1/ vdu nor the uncertainty in deter- 
mining p c . Moreover, logarithmic corrections to scaling 
should be considered to obtain accurate values of the ex- 
ponents. 



IV. CONCLUSION 

The first important fact is to notice that critical expo- 
nents and scaling relations obtained by simulations are 
globally compatible with the expected theoretical values 
calculated for bond percolation. This gives confidence in 
the extraction of the (unknown) site percolation thresh- 
old. This also confirms universality between site and 
bond percolation for this model. 

However, the main result of this paper is the value 
p c = 0.7360(5) for site percolation on planar <& 3 random 
graphs. It is greater but not very far from the thresh- 
old value for site percolation on the honeycomb lattice 
p c (honeycomb) = 0.6962..., which is the simplest regular 
trivalent lattice. On one hand, this means that, to some 
extent, planar $ 3 random graphs and honeycomb lattices 
look alike: more precisely , they locally look alike. In con- 
trast, trivalent Bethe lattices are neither locally equiv- 
alent to honeycomb nor to planar <I> 3 random graphs, 
so that their percolation thresholds are very different, 
p c (3 — Bethe lattice ) = 1/2. On the other hand, as p c 
is greater for planar <& 3 random graphs than honeycomb 
lattices, percolation is easier on a pure hexagonal lattice 
than on the planar $ 3 random graphs. The reason is that 
on these latter graphs, there are regions called baby uni- 
verses (B.U.) connected to the rest of the graph by very 
small boundaries called necks Moreover, B.U. can 
grow on other B.U., giving a fractal (self-similar) struc- 
ture to the graph. So, for a given occupation probability 
p, the probability that a given B.U. belongs to a giant 
connected cluster is proportional to the probability that 
at least one vertex of its boundaries is occupied. This is 
small compared with the probability that, on a pure hon- 
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eycomb lattice, a given region is a part of a giant cluster. 
This fractal structure of B.U. is also the main feature that 
makes honeycomb lattice and planar <I> 3 random graphs 
globally different at long distance, so that their critical 
exponent are different. 

It would be interesting to study in more details the con- 
nections between baby universes and percolation tran- 
sition. In particular, a non-uniform occupation proba- 
bility, depending for instance on the local curvature or 
on the structure of B.U., could shed light on this prob- 
lem. The role of B.U. could also be studied by real-space 
renormalization group analysis. As mentioned above, 
planar <£> 3 random graphs have a hierarchical structure 
that makes them look like trees of baby universes. It 
is possible to use this self-similarity of planar $ 3 ran- 



dom graphs with respect to B.U. to perform a real-space 
renormalization group transformation [l5[ by replacing 
each baby universe of last generation (i.e. a B.U. with no 
further B.U. growing on it) by one supersite Then, 
if the corresponding last generation B.U. (including its 
boundary) contained a spanning cluster, the supersite is 
occupied. This defines an occupation probability p' for 
the supersite as a (complicated) function of p, the occu- 
pation probability of the original graph. 

It should also be noticed that the value of p c found here 
is comparable with high values found on Archimedean 
lattices [lfj]. It would be interesting to understand if pla- 
nar $ 3 random graphs share common local characteristics 
with Archimedean lattices. 
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